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Abstract 

We present the solution of a linear Restricted Solid-on-Solid (RSOS) model in a field. 
Aside from the origins of this model in the context of describing the phase boundary in a 
magnet, interest also comes from more recent work on the steady state of non-equilibrium 
models of molecular motors. 

While similar to a previously solved (non-restricted) SOS model in its physical behaviour, 
mathematically the solution is more complex. Involving basic hypergeometric functions 
3</>2, it introduces a new form of solution to the lexicon of directed lattice path generating 
functions. 

1 Introduction 

The SOS model arose from the consideration of the boundary between oppositely magnetised 
phases in the Ising model [11 at low temperatures and is now considered to be useful for describing 
the salient features of a wide variety of interfacial phenomena [2j |3j HJ |5j |6] . The configurations 
involved in the linear (1+1) -dimensional case, modelling the interface in a two-dimensional 
magnet, have also been used to model the backbone of a polymer in solution. The critical 
phenomena associated with this model describe wetting transitions of the interface with a wall 
[2]. For the SOS model the phase diagram contains a wetting transition at finite temperature 
T w for zero field and complete wetting occurs taking the limit H —* for T >T W [7] . 

The linear SOS model with magnetic field and wall interaction was solved in |7J. The 
restricted SOS model is a variant of the SOS model where the interface takes on a restricted 
subset of configurations as opposed to the full SOS model. Effectively it suppresses very large 
local fluctuations of the interface. This model has been considered with wall interaction, but 
not yet with magnetic field. This is partially because, with wall interaction only, the SOS 
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and RSOS models are in the same universality class, as demonstrated by their exact solutions 
IS El HOI ED H21 E3]. Recently, it has been suggested that the RSOS model in a field describes 
the steady state of a non-equilibrium model of a molecular motor [T3]. This observation has 
motivated us to consider the RSOS model in a field. In doing so, we have discovered a novel 
form of generating function for a directed lattice path problem, and a new method of solution 
for such problems. 

The RSOS model we analyse can be described as follows. Consider a two-dimensional square 
lattice in a half plane. For each column i of the surface a segment of the interface is placed on 
the horizontal link at height Ti > and successive segments are joined by vertical segments to 
form a partially directed interface with no overhangs. The configurations are given the energy 

-/3E = -Kj2\n-ri-i\-Hj2ri + bJ2 5 n,o ■ (1-1) 

i i i 

There are two basic variants of this model which have been discussed in the literature. If there 
are no restrictions placed on the differences of successive heights rj, the model is called the 
(unrestricted) SOS model, analysed in [7]. 

On the other hand, constraining the height differences to be bounded by 

| Ari| = h-r;_i| < 1 (1.2) 

gives the restricted SOS (RSOS) model. This has previously only been considered in the case 
of zero field H [8j [9], for several types of external potential [15]. Both variants have been 
considered, utilising a different thermodynamic ensemble, as models for polymers in solution, 
since the finite configurations are partially directed self-avoiding walks [8j [16] . In [T7J a RSOS 
model with H = but a rigidity term dependent on \An — Arj_i| has been considered as a 
model of semi-flexible polymers such as DNA. 



2 The RSOS generating function 

We discuss the RSOS model in terms of lattice paths. An RSOS path is a partially directed 
self-avoiding path with no steps into the negative x-direction and no successive vertical steps. 
To be precise, an RSOS path of width N with heights ro to r^ has horizontal steps at heights 
ri, . . . , rjVj an d vertical steps between heights rj_i and r, for i = 1, . . . , N, but no horizontal step 
associated with ro. This means that an RSOS path starts at height ro with either a horizontal 
step (if r\ = ro) or vertical step (if r\ ^ ro), but must end at height r/v with a horizontal step. 
Figure [T] shows an example. 

The partition function for the RSOS paths of width N with ends fixed at heights ro > and 
fN > 0, respectively, is given by 

„, , I exp(-/3£(r ;ri)) , \r - n\ < 1 

Zx{ro;n) = < (2.1) 
I , |r - ri| > 1 , 
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Figure 1: A typical RSOS configuration beginning on the surface and finishing on the surface 
with a horizontal step: each horizontal step is assigned a weight x, each vertical step a weight 
y, each unit of area a weight q, and each step that touches the surface an additional weight 
k. The width of the configuration shown here is N = 9, and the heights are ro = 0, t\ = 1, 
V2 = ?~3 = r4 = 2, rs = 1, r§ = 0, r-j = rs = 1, and rg = 0. The weight of this configuration 
equals x 9 y e q 10 K 2 . 



and 

Z N (r ;r N )= ^ exp(-/3£(r ; n, . . . , r N )) , N = 2, 3, . . . , (2.2) 

ri,...,rjv_i>0 
In— »"i-l|<l 

where 

N N 

-(3E(ro;n,...,r N ) = -Kj^ \n - n~i\ - + bj^tnfi ■ (2.3) 

i=l i=l i=l 

Here, we shall consider paths with both ends attached to the surface, i.e. we shall focus on the 
partition function 

Z N = Z N (0;0). (2.4) 

We define 

y = exp(-K) , q = exp(— H) , and k = exp(6) , (2.5) 

so y is a temperature-like, q a magnetic field-like and n a binding energy-like variable, and 
write 

Z N = Z N (y,q,n) . (2.6) 

The free energy is then 

- (3f(y, Q, «) = Jim log Z N (y, q, k) . (2.7) 

AT— >oo JS 

Define the generalised (grand canonical) partition function, or simply generating function, as 

oo 

G(x,y,q,K) = l+^2x N Z N (y,q,K) . (2.8) 
JV=1 

Thus, the radius of convergence x c (y, q, k) of G(x, y, q, k) with respect to the series expansion in 
x can be identified as exp(/3/(y, q, k)), hence 

f(y,Q,K) = kT log x c (y,q,K) . (2.9) 
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It is convenient to consider G as a combinatorial generating function for RSOS paths, where 
x, y, q and k are counting variables for appropriate properties of those paths. Interpreted in 
such a way, x and y are the weights of horizontal and vertical steps, respectively, q is the weight 
for each unit of area enclosed by the RSOS path and the x-axis, and n is is an additional weight 

for each step that touches the surface. For example, the weight of the configuration in Figure [T] 

q fi in 2 
is x y q k . 

We find easily the first few terms of G as a series expansion in x, 

G{x, y,q,K) = l + KX + (k 2 + Ky 2 q)x 2 + . . . , (2.10) 

where the constant term corresponds to a zero-step path starting and ending at height zero with 
weight one. 

3 Functional equation and exact solution 

The key to the solution is a combinatorial decomposition of RSOS paths which leads to a 
functional equation for the generating function G. This decomposition is done with respect to 
the left-most horizontal step touching the surface at height zero, and is shown diagrammatically 
in Figure |2j 




G(x,k) = 1 + kxG(x,k) + Kx 2 y 2 qG(qx, l)G(x, k) 

Figure 2: The diagrammatic form of the functional equations for RSOS paths, indicating the 
combinatorial decomposition of RSOS paths. 

We distinguish three cases: 

(a) The RSOS path has zero width, and there is no horizontal step at height zero. The 
contribution to the generating function is 1. 

(b) The RSOS path starts with a horizontal step, which therefore is at height zero. The 
rest of this path is again a RSOS path. The contribution to the generating function is 
kxG(x, y, q, k). 

(c) The RSOS path starts with a vertical step. Then there will be a left-most horizontal 
step at height zero, and removing this step cuts the path into two pieces. The left path 
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starts with a vertical and horizontal step, followed by an RSOS path starting and ending 
at height one and not touching the surface, followed by a vertical step to height zero. 
The right path is again a RSOS path. The contribution to the generating function is 
yxqG(qx, y, q, l)ynxG(x, y, q, k). 

Put together, this decomposition leads to a functional equation for the generating function 

G(x, y, q, k) = 1 + kxG(x, y, q, «) + KX 2 y 2 qG(qx, y, q, l)G(x, y, q, k) . (3.1) 
For q = 1, the solution is a simple algebraic function 



G(x,y,l,K) 



1 — KX 



1 — x — y/ (1 — x) 2 — 4x 2 y 2 



(3.2) 



and for general values of q, a formal iteration of (3.1 ) leads to a continued fraction expansion 

1 



G(x,y,q,K) 



1 — KX 



2 1 

nqx y 



(3.3) 



1 — qx 



q 3 x 2 y 2 



1 



q 5 x 2 y 2 



1 — q 3 x 



fx 2 y 2 



1 — q 4 x — . . . 

However, there is a non-trivial method to solve the functional equation for G in terms of power 
series. Our main result is an expression for G involving g-series. In the next section, we shall 
derive the following expression: 



G(x,y,q, 1) 



a(Xqx;q)oo(p 




+ ((1 


- X)qx; q)oo4 


'[(1 




a(Ax;g)oo0 


^x,qj^;q,q 


+ ((1 


- X)x; q)oo(j) 


[a- 





with 



and 



a 



[0,q^;q,q] - (1 - X)cf> [0, g±^; q, q 



Xcf) 



Q,qj^x;q,q 2 



y 2 = A(i - A) 



0, qj^x]q,q 



Here <f> is given in terms of the basic hypergeometric function 3^2 as 



, . .'0,0,0 
[s.l:q. :J - ;-;() 2 j s t ;q,z 



E 

ra=0 



(s;q) n (t;q)n(q;q)n 



and 



n-l 



(t;g)n = n( 1 -^) 



k=0 



is the standard (/-product. From Equation (3.1) we have the full solution as 

G{x,y,q,K) 



1 — kx — KX 2 y 2 qG(qx, y, q, 1) 



(3.4) 



(3.5) 



(3.6) 



(3.7) 



(3i 



(3.9) 
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4 Solving the functional equation 



Using a linearisation Ansatz [18] , standard for a g-deformed algebraic equation such as Equation 
(3.1), we substitute 

G(x,y, g ,l)= R{x ^ q) (4.1) 



into Equation (3.1) with k = 1, and find that H(x,y,q) must satisfy the linear g-functional 
equation 

H (x, y, q) + {x- l)H(qx, y, q) + qx 2 y 2 H(q 2 x, y, q) = . (4.2) 
We then try to solve this linear functional equation using a series in x, 

oo 

H(x,y,q) =Y,c n (y,q)x n . (4.3) 

n=0 

This unfortunately leads to a non-trivial three-term recurrence for the coefficients 

(1 - q n+2 )c n+2 (y, q) + q n+1 c n+1 (y, q) + y 2 q 2n+1 c n (y, q) = (4.4) 

for n > with initial condition co(y, q) + (1 — q)ci(y, q) = 0. This is different from the two-term 
recurrences obtained when considering the models in [18] which can be solved by direct iteration. 

Inspired by the structure of basic hypergeometric functions, which we know form the basis 
for a solution in the SOS model [7], we transform the coefficients as 

/_ ^\riqn{n- 1)/2 

c n (y,q) = -r—, d n (y,q) . (4.5) 

This leads to a recurrence 

d n +2{y, q) - d n+ i(y, q) + y 2 d n (y, q) = y 2 q n+1 d n (y, q) (4.6) 

for n > with initial condition d\{y,q) = do(y,q). While this is still a three-term recurrence, 
only one of the terms has a non-constant coefficient. 



The left-hand side of Equation (4.6) is a homogeneous difference equation with constant 
coefficients and characteristic polynomial 

P(A) = A 2 -A + y 2 . (4.7) 



If the right-hand side of Equation (4.6 ) were zero, the solution would be given as d n = AX^+BX 1 ^ 
where Aj are the roots of P{X) = 0. To solve this recurrence, we use the Ansatz |8j 

oo 

d n (y,q) = X n Y,e m (y,q,X)q nm ■ (4.8) 

m=0 



Inserting this Ansatz into Equation (4.6) we find 



P(A)e (y, q,X) + J2 ? m " [P^q m )e m (y, q, A) - y 2 qe m ^(y, q, A)] = . (4.9) 



m=l 
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Necessarily P(X) = 0, and normalising by letting eo(y, q, A) = 1, we find by iteration 

q m 



e m {y,q,\) 



{X 2 q/y 2 ;q) m (q;q) m ' 



(4.10) 



The full solution to the recurrence equation (4.6) is a linear combination of (4.8) over both 
values of A satisfying P(X) = 0. If P(X) = 0, then also P(y 2 /X) = 0, and we can write 



d n (y, q) = AX n J2 em{v, q, X)q nm + B(y 2 /X) n £ e m (y, q, y 2 /X)q r ' 



(4.11) 



m=0 



m=Q 



where A is an arbitrarily chosen solution of A(l — A) = y 2 . Using the initial condition do(y, q) = 
di(y, q) we can solve for the ratio a = A/B. We can somewhat simplify our expressions by noting 
that the dependence between A and y implies that we can replace y 2 by A(l — A). Noticing that 
the functions involved are 291 basic hypergeometric series, we have 



201 ;q,qj-(l- A) 2 0i ( g ^A ; q, q 2 ^ 



A 20i ( „°^_;q,q 2 

q 1-A 



29l 



A 

0,0 

«1-A 



(4.12) 



where 



29l 



0,0 



( s . 



n=0 



{s;q)n(q;q)n 



3<P2 



0,0,0 

0,8 ' 



(4.13) 



Substituting the expression for d n given in (4.11) back into H using Equations (4.5) and (4.3), 
we find that 



H(x,y, 



-a 



E 

m,n=0 
00 

E 



1 n(n—l)/2/ 



-Xx) n q 



n n nm+m 



(q\/ (1-A); q) m (q; q) m (q; q) n 
qn(n—i)/2/n _ \^x) n q nm+ ' n '^ 
(q(l - A)/A; q) m (q; q) m {q; q), 



(4.14) 



m,n=0 

The summation over n can be done explicitly using Euler's formula [T9 

^ra(n— l)/2/ 



(*;?)oo = Yj 



n=0 



(q;q)r 



(4.15) 



We find that 

00 

H(x,y,q) =aJ2 



(Xxq^q)^ 



m=0 



(qX/(l - A); q) m (q; q) 



^ (a(l 



{{1-X)xq m ;q) 00 q r 



— (?(1-A)/A;g) m (g;g), 



(4.16) 



which after pulling out a (/-product factor from each sum can be identified with 

0,0,0 



H(x, y, q) =a (Ax; <?)oo3</>2 



Xx,qX/(l - A) 



;q,q 



(4.17) 



+ ((1 - X)x;q) oo3 <j)2 



0,0,0 
(1-A)x,g(l-A)/A 



q,q 



Substituting into the linearisation Ansatz (4.1) gives us the solution written in Equation (3.4). 
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5 Comparison of SOS and RSOS model solutions 



The unrestricted SOS model was solved in [7] by a different technique, namely the Temper ley 
method. For the sake of comparison it is worthwhile reproducing the solution of the SOS 
model via the functional equation technique presented above for the restricted SOS model. This 
highlights a fundamental difference in the difficulty of solving the two models. 



In analogy to the functional equation (3.1 ), the SOS model generating function S(x, y, q, k) 
satisfies 



S(x,y,q,n) = 1 + k x(l - y 



q 



S(x,y,q,K) + K—S(qx,y,q, l)S(x,y,q,n) 



(5.1) 



where the variables have identical meaning. 

Setting k = 1 and substituting the linearisation Ansatz S(x,y,q,l) = T(qx,y,q)/T(x,y,q) 



analogous to the one used above in Equation (4.1) yields 

;2- 



T(x,y,q)+(x(l-y z )-l 



y 



T{qx, y, q) H T(q x,y, q) = . 



(5.2) 



While superficially this may seem very similar to the linear functional equation (4.2) for H, it 
is important to recognize that the method of solution is to substitute a series in the variable x 



into the functional equation. Equation (5.2) above contains only linear factors in x, while the 



functional equation (4.2) contains factors that are quadratic in x, This has the effect that the 
resulting recurrence for the coefficents of the series expansion of T in x is a two-term recurrence, 
while we found a three-term recurrence in the case of the RSOS model. The two-term recurrence 



leads to an immediate solution by iteration, while the three-term recurrence (4.4) requires the 
extra work detailed in the previous section. One readily finds 

(_l)« g n(«-l)/2[ x (l _ y 2j]n 

n=0 



T(x,y,q) = ^2 



(y 2 ;q)n(q;q) n 



191 



j,x(l - y 2 ) 



(5.3) 



which leads to an expression equivalent to Equation (45) in [7J. The resulting expression for 
T should be compared with the much more complicated expression for H given in Equations 



(4.17) and (4.12). 



6 Conclusion 

In this paper we have presented a solution to the linear RSOS model in a field. We have 
expressed the solution in terms of basic hypergeometric functions at values of their arguments 
which are not powers of the counting variable in a combinatorial problem. We know of one other 
case where this occurs in a different manner [20 . 

We note that there is a bijection between RSOS paths and Motzkin paths. Hence it would 
be interesting to consider other lattice path problems with a similar structure such as fc-coloured 
Motzkin paths [21] . 

It will also be interesting to see if this solution will produce insights into non-equilibrium 
models of molecular motors. 
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